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Abstract 

In this paper, we study the impact of frequency-selectivity on the gap between the required per-user 
transmit power in the residual CFO scenario (i.e. after CFO estimation/compensation at the base-station 
(BS) from lb]) and that in the ideaFzero CFO scenario, for a fixed per-user information rate, in single¬ 
carrier massive MU-MIMO uplink systems with the TR-MRC receiver. Information theoretic analysis 
reveals that this gap decreases with increasing frequency-selectivity of the channel. Also, in the residual 
CFO scenario, an 0{'/M) array gain is still achievable (M is the number of BS antennas) in frequency- 
selective channels with imperfect channel estimates. 

Index Terms 

Massive MIMO, carrier frequency offset (CFO), frequency selective channel. Time-reversal maxi¬ 
mum ratio combining (TR-MRC), single-carrier. 

I. Introduction 

Massive multi-user (MU) multiple-input multiple-output (MIMO) system/large seale antenna 
system (USAS) has been identified as one of the key next generation wireless teehnologies, due to 
its eharaeteristie ability to provide huge inerease in energy and speetral effieieney with inereasing 
number of base-station (BS) antennas ([IJ. In a massive MU-MIMO system, the USAS present 
in the BS serves an uneonventionally large number of single-antenna user terminals (UTs) in 

Sudarshan Mukherjee and Saif Khan Mohammed are with the Department of Electrical Engineering, Indian Institute of 
Technology (I.I.T.) Delhi, India. Saif Khan Mohammed is also associated with Bharti School of Telecommunication Technology 
and Management (BSTTM), I.I.T. Delhi. Email: saifkmohammed@gmail.com. This work is supported by EMR funding from 
the Science and Engineering Research Board (SERB), Department of Science and Technology (DST), Government of India. 


2 


the same time-frequeney resouree O. It has been shown that for a given number of UTs in a 
coherent massive MU-MIMO system, with imperfeet ehannel estimates, the required per-user 
transmit power (to aehieve fixed per-user information rates) ean be reduced as with increasing 
M (i.e. 0{^/M) array gaiio, where M is the number of BS antennas 0. 

The result discussed above assumes perfect frequency synchronization between the BS and the 
UTs for coherent multi-user communication. In practice, acquisition of carrier frequency offsets 
(CFOs) between the carrier frequency of signals received from different UTs and the BS oscillator 
is a challenging task in massive MIMO systems (because of unconventionally large number of 
UTs). It has been observed that the existing optimal/near-optimal CFO estimation techniques 
for small scale MIMO systems are not amenable to practical implementation in massive MIMO 
systems due to prohibitive increase in their complexity with increasing number of UTs and also 
with increasing number of BS antennas flSl, 0. In jS], the authors consider an approximation 
to the joint maximum likelihood (ML) estimator for CFO estimation in frequency-flat massive 
MU-MIMO systems and analyze the information rate performance. However the CFO estimator 
presented in [|5l requires multi-dimensional grid search and hence has high complexity with 
increasing number of UTs. Also the information rate analysis in dH does not consider the 
imperfect CSI scenario. Recently in [[61, a low-complexity near-optimal CFO estimation technique 
has been proposed for massive MU-MIMO systems. Using this technique in Q, the information 
rate performance in frequency-flat Rayleigh fading channel has also been studied. However for 
frequency-selective single-carriecl massive MU-MIMO uplink systems, it is not known as to how 
the information theoretic performance gap between the ideal/zero CFO scenario and the residual 
CFO scenario (i.e. with CFO estimation/compensation from 0) would vary with increasing 
frequency-selectivity. In this paper, we study this performance gap and show the interesting 
result that it decreases with increasing frequency-selectivity. 

The novel contributions of our work presented in this paper are as follows: (i) we have derived 

'Under the average power constraint, for a fixed M, fixed number of UTs and a fixed desired per-user information rate, the 

per-user average transmit power decreases with increasing M m 

^It has been shown that a single-carrier massive MU-MIMO system can achieve performance comparable to their OFDM 


counterparts i). 
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Fig. 1 The communication strategy: CFO Estimation and Compensation Strategies and Data Commipiication. 

a closed-form expression for per-user information rates of eaeh UT with the TR-MRCO reeeiver 
at the BS in the imperfeet CSI seenario, in the presenee of residual CFO errors (after CFO 
eompensation using the CFO estimator in [[6l|); (ii) analysis of the information rate expression 
reveals that an 0{'/M) array gain is still aehievable in frequency-selective channel with CFO 
estimation/eompensation, i.e., no loss in array gain when eompared to the ideal/zero CFO 
seenario; (iii) further analysis of the information rate expression reveals the interesting result that 
for a fixed per-user information rate, the gap in the required per-user transmit power between 
the residual CFO seenario and the ideal/zero CFO seenario decreases with increasing frequency- 
selectivity (i.e. number of ehannel memory taps L) of the ehannel. For instanee with M = 160 
BS antennas, iC = 10 UTs and a per-user rate of 3 bits per ehannel use (bpeu), this gap in per¬ 
user transmit power is approximately 4.22 dB when L = 1 and is approximately 0.07 dB when 
L = 20. [Notations: E denotes the expeetation operator. (.)* denotes the eomplex eonjugate 
operator.] 

II. System Model & CFO Estimation 


We eonsider a single-eell single-earrier massive MU-MIMO BS with M antennas serving 
K single antenna UTs in the same time-frequeney resouree. The baseband frequeney seleetive 
ehannel is modelled as a diserete time finite impulse response (FIR) filter with L ehannel taps. 
The ehannel gain eoeffieient from the kf' UT to the m* BS antenna at the ehannel tap is 
given by hmk[l] ~ where {(Jhki > 0} models the power delay profile (PDF) of 

the ehannel for the fc* UT {k = 1, 2,... , J’f; I = 0,1,. .., L — 1 and m = 1, 2, ..., M). Sinee 
a massive MU-MIMO system is expeeted to operate in time-division duplexed (TDD) mode, 
eaeh eoherenee interval {Nf. ehannel uses) is split into an uplink (UL) slot {N^ ehannel uses) 

^TR-MRC (Time-reversal maximum ratio combining) receiver is a well-known low-complexity single-carrier multi-user 
detector in frequency-selective massive MIMO channels E). 
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and a downlink (DL) slot (N^ — Nu channel uses). From the eoherent eommunieation strategy 
depleted in Fig. [H it is observed that CFO estimation is performed in a speeial UL slot prior to 
the conventional UL eommunieation (channel estimation/UL data transmission). 

For the CFO estimation phase, the UTs transmit pilots having average power for N < 
Nu ehannel uses. The pilot sequenee is assumed to be divided into B = \N/KL] bloeks 
of KL ehannel uses eaeh. The UT transmits an impulse of amplitude y/KLpa at time 
t = [b — 1)KL + {k — 1)L + I and zero elsewhere (b = 1, 2,..., i?; I = 0,1,..., L — 1 
and k = 1,2,..., K). Using bloek-wise eorrelation of the reeeived pilots, CFO estimation is 
performed at the BS, using the low-eomplexity CFO estimator proposed in 10. Let cuk be the 
CFO of the /c* UT and cD^ be the estimate of cUk- It can be shown from the eentral limit 
theorem (CUT), that as M —)■ cx), the error in CFO estimation is asymptotieally Gaussian, i.e., 
(cDfe — Uk) ~ ^ ~ is the mean squared error (MSB) of CFO 

estimation for the /c* UT and is given by 0 

Ik V B-l 2K^k 


^Ik = 


wherein jk = ^ J2f=o ^Iki i® the reeeived SNR from the /c* is the AWGN power at the BS) 


M{N - KL){KLfGl 


2 ’ 


( 1 ) 


and Gk = Em=i \hmk[l]?/{M (^Iki)- Note that with M cx, from the strong law 

of large numbers, for independent hmk[l], lim Gk = 1- 

M^OO 

Remark 1.' (Proposition 2 from From the expression of in dU) it is elear that with 
increasing M ^ oo, fixed N, K and L, the required received SNR 7 ^ (to achieve a fixed desired 

MSB of CBO estimation) deereases as -4^, i.e., with 7 ^ cx we have lim = eonstant. 

^ Cm’ Cm’ m^oo 

□ 


^L-l 




III. Uplink Data Communication 


After the CBO estimation phase, the eonventional UB data eommunieation starts at f = 0 (see 
Big. [T]). Note that the CBO eompensation is performed at the BS prior to UB channel estimation 
and also prior to UB receiver proeessing. Bor UB data eommunieation, the first KL ehannel 
uses are dedieated for UB pilot transmission, which is then followed by (L — 1) ehannel uses 

^Note that the low-complexity CFO estimator in (6) is well defined if and only if \ujkKL\ <C tt. For most practical massive 
MIMO systems, this condition would hold true (6). 
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M L-1 L-1 


ISl4t] = Y1 ^*rnk[l] hmk[l'] Xk[t -/' + /] (2) 

m=l 1=0 l'=0,l'^l 

M L-1 K L-1 

MUIfc[t] = y/p^ EE E E Knk [/] h^,[l'] x,[t -I' + l] e^'((2,-s.)(t+0-A<..(t-(g-i)L+(«-z'))) (3) 

m=l 1=0 q=l,q^k l'=0 


E 


ESk[t]W:[t] =E[A4 


E 


Xk[t]\ Al[t]xl[t] - E[Ak[t]]xl 


+ E 


ISl4t]+MUI, 


= 0, since X}^ [t] and [t] are independent 


= 0, since all Xk[t] are i.i.d. 


E 


Xfc[t]EN*[t] 


= 0 . 


(4) 


= 0 (since alln^fc[£]are i.i.d., 
zero mean and independent 
of Xklt]) 


of pre-amble sequence. The UL data transmission occurs in 
followed by (L — 1) channel uses of post-amble transmission 


le next channel uses and is 


A. Channel Estimation 


In the channel estimation phase, the UTs transmit sequentially in time, i.e., the kA UT transmits 
an impulse of amplitude yJKLpa at time t = {k — 1)L and zero elsewhere. The received pilot at 
the m* BS antenna at f = {k — l)L+l is given by rm[{k — l)L+l] = y/KLpa hmk+ 
Wm[{k — 1)L + T\, where m = 1, 2,..., M; I = 0,1,..., L —1 and k = 1,2,..., K. Here Wm[ik — 
1)L -I- /] ~ CAf{0, a^) is the circular symmetric AWGN. We first perform CFO compensation 
for the fc* UT by multiplying rm[(fc — 1)L -f /] with which is then followed 

by computing the maximum likelihood estimate of the channel gain coefficient, i.e^ hmk[l] = 
-j^rm[{k - l)L + + -^^^nmk[{k - 1)L + 1], where^ hmk[l] = 

~ CJ\f{0, is the effective channel gain coefficient and nmk[{k — l)L + 
1] = Wm[{k — 1)L + ~ CJ\f{0,a'^). Here Auk = Qk — ^k is the residual CFO 


error. 

^The symbols transmitted in the pre- and post-amble sequences are independent and identically distributed (i.i.d.) and are 
assumed to have the same distribution as the information symbols (see section III-B) to ensure the correctness of the achievable 
information rate expression. 

®Both hmk[l] and Wm[{k — 1)L + 1] have uniform phase distribution (i.e. circular symmetric) and are independent of each 
other. Clearly, rotating these random variables by fixed angles (for a given realization of CFOs and its estimates) would not affect 
the distribution of their phases and they will remain independent. Therefore the distribution of hmk[l] and nmk[{k — 1)L + 1] 
would be same as that of hmk[l] and Wm[{k — 1)L -f 1] respectively. 
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B. Uplink Receiver Processing 


After channel estimation and preamble transmission, the UL data transmission begins at time 
t = KL+L — 1. Let Xk[t] ~ CAf{0, 1) be the i.i.d. information symbol transmitted by the k-th UT 
at the f-th channel use and pu be the average per-user transmit power. The received signal at the m- 


th BS antenna at time t is therefore given by rm[t] = ^/IhJ2q=l J2i=o hmqmXq[t — l]e^‘^'>^ + Wm[t], 
where t = KL + L — 1,..., KL + L + Nqj — 2(= — L). To detect Xk[t], we first perform 

CFO compensation for the UT on the received signal, followed by TR-MRC processing |0. 
Output of the TR-MRC receiver at time t is given by 

M L-1 

Mt] = Y, + ^ Ak[t] Xk[t] + ISIfc[t] + + ENfc[f], (5) 

^ ^ CFO Compensation ^ 

where step (a) follows from the expressions of h^kil] (see Section III-A) and r^lt] (see Section 

M L-l _ 

III-B) and Ak[t] = ^ ^ Here the terms ISIfc[t] (inter-symbol 

m=l Z=0 

interference), MUIfc[f] (multi-user interference) are given by Q-® at the top of the last page 


and ENfc[f] = Xk[t] — Ak[t]xk[t] — ISIfc[f] — MUIfc[f] (xk[t] is given by the first line of ®). From 

( 6 ) 


®, we therefore have 

Xk[t] = IE[Afc[f]] Xk[t] + SIFfc[f] -f ISIfc[f] -f MUIfc[f] -f ENfc[t], 


= ESfc[t] 


= Wk[t] 


where SIEfc[f] = (Afc[f] — E[Afc[f:]])xfc[f] is the time-varying self-interference component and 
ESk[t] is the effective signal componentO Note that the statistics of both ESk[t] and Wk[t] are 
functions of t. However for a given t, the realization of Wk[t] is i.i.d. across multiple UE data 
transmission blocks (i.e. coherence intervals). Therefore for each t, we have a SISO (single-input 
single-output) channel in ®, when viewed across multiple coherence intervals. Thus for Njj 
channel uses, we have Nn SISO channels with distinct channel statistics. We therefore have Njj 


different channel codes, one for each of these Nd channels. ' 
use of every coherence interval is jointly decoded at the BSH 


le data received in the channel 


C. Achievable Information Rate 

Since Xk\f\ and nmk\p] are all independent and zero mean, it can be shown that ESfc[f], SIEfc[t], 
ISIfc[f], MUIfc[f] and ENfc[f] are all zero mean and uncorrelated with one another (see @ at the 

^In E[.] is taken across multiple channel realizations and also across multiple CFO estimation phases. 

*This coding strategy has also been used in Q, ©. 
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top of the page). Since the overall noise and interference term Wk\t] and ESfc[f] are uncorrelated, 
a lower bound on the information rate for the effective channel in ® can be obtained by 
considering the worst case uncorrelated additive noise (in terms of mutual information). With 
Gaussian information symbols Xk\t], this worst case uncorrelated noise is also Gaussian with the 
variance E[|W4t]|2] = E[|SIF4f]|2] +E[|ISl4f]|2] +E[|MUl4t]|2] +E[|EN4f]|2] [[TOl. Therefore 
an achievable lower bound on the information rate for the channel is given by I{xk[t]] Xfe[f]) > 
log 2 (l + SINRfe[f]), where SINRfc[f] = E[|ESfc[f]p]/E[|H4[f]p]. Therefore the information rate 
for the fc* UT is given by 1^ = + SINRfc[f]). Using the expressions of 

ESfc[f], SIFfc[f], ISIfc[t], MUIfc[f] and ENfc[f] from Q,®, ® and ®, and Accfc = {uk — ojk) ~ 
AA(0,(T^^), the variances of each term are computed and summarized in Table HI Using Table HI 
SINRfc[f] = E[|ESfc[t]n/E[|Wfc[f]|2] is given by 


SINRfc[f] = 




[1 _ I 1 I Cl I C2 ’ 

^ J MK-y? ^ M'vt. ^ M 


(7) 


where ci = 1 + 


MK^l 
A 


and C 2 = j- Ell Oq- Here Ok = Yto ^hki- 


M'Yk 

^2 


Remark 2: (Array Gain) In the following, for a fixed desired information rate, and therefore 
fixed SINRfc[f], we examine the rate of decrease in the required jk with increasing M. On the 
RHS of © we note that the numerator and the first term in the denominator depend on the 
received SNR yk only through the MSE of CEO estimation, . Erom Remark [T] we know that 


TABLE I List of Variance of all components of Wfc[d- 


Component 

Variance 

ESfc[d 

A/r2 2 U (t-(fc-l)Z,)2 

M Pul E ^hki 1 e “C ' 

^ i=0 ' 

SIFfe[t] 

L-l 

+ Mpu E ^hkl 
i=0 

ISIfc[d 

r / -L-l , 2 L-l , 

^ (^ hkl ) ~ ^ hkl \ 

^^ 1=0 ^ (=0 ^ 

MUL[d 

/L-l V / K L-l V 

( 51 51 ^hqVj 

^ 1=0 ' ^ q = l , qj^k L = 0 ' 

EN;,[d 

^ f E +^ + Ma ^ (llalY 

\q = l 1 = 0 J ^ \l = 0 J 













if 7 fc oc then as M —)■ oo, converges to a constant, i.e., the numerator and the first term 
in the denominator of ([7]) converge to constant values. Further as M —)■ oo, the last term in the 
denominator of (Q, i.e, ^ vanishes. The rest two terms, ( 77 ^ + ]^) however depend on 
both M and 7 ^. Since with decreasing 7 ^, the term 77^-7 would eventually dominate the other 
term we must therefore decrease 7 ^ as so that the SINRfc[f] converges to a constant as 
M —>■ 00 (This shows that with every doubling in M, 7 ^ decreases roughly by 1.5 dB for a fixed 
per-user rate when M —)■ 00 (see the change in 7 ^ from M = 320 to M = 640 in Table HH). □ 


Lemma 1: Consider \uJkKL\ <C vr and lim M 7 ? = constant > 0. With fixed K, N and a 

M^OO 

fixed desired information rate of the channel code Rk[t] = lim log 2 (l + SlNRfc[f]) Rq k[t] 
(where Ro,k[t] = log 2 (l + ak,t), ak,t = 2 (jv-^^L)(XL)^ )’ asymptotic (i.e. M ^ 00 ) gap 
between the required 7 ^ in the residual CFO scenario (i.e. after CFO estimation/compensation) 
and that in the ideal/zero CFO scenario (jk o)> i-e-> lim — decreases with increasing L, 

’ M^oo 

provided L < 


Proof: Since lim M 7 ? = constant > 0 and lim Gk = 1, from (fTl) we have lim (t — 

M^oo M^oo M^-oo 

lim —j(t—(fc—1)L)^ , 

{k - 1)LY = '^^ 2 (iv-^L)(KL)^ -= Oik,t0k, where 9k = lim Using this limit in O, we 


have Rk[f\ = lim log 2 (l + SlNR^lfl) = log 2 (l + e (1 — + Ok)). Let 6k = 9' be the 

M^oo 


unique solution to this equation for a fixed i.e., 

(1 + 4)(1 - 2 ”^'=^) = 

Therefore (1 + 9'){1 — 2“'^'=^) = < 1 9' < 


( 8 ) 


1 






2 -^ 0 ,A; 


1 

^k,t 


ak,t9' < 1 


o ^k.t^k 


1 — Oik,t9k. Substituting this in ([ 8 ]) with 9k = 9', we have = 


(^k,t + (1 + ^) ~ (1 + C.' Rk\P\ S> Ro,k[t] > log 2 (l + + otk^t)))- 

Similarly for the zero CFO scenario, using = 0 in (|7]) we have 9q = lim 


2 «aW-1 - 


^The system parameters for data in Table [II] is given in the first paragraph of Section IV. 


TABLE II Minimum required 7 ^ (in dB) for fixed information rate /fc = 1 bpcu (k = 1) with increasing M, K = 10 
and L = 10 


M = 40 

M = 80 

M = 160 

M = 320 

M = 640 

-9.5266 

-12.2927 

-14.5049 

-16.4462 

-18.2341 






















From the expressions of 6o and 6', the asymptotie SNR gap is given by 


lim = lim 


M- 


°° 7a:,o m^ooV 1/MK'y^ 


mKii, ^— 

= v -^ = V1 


( 9 ) 


Since ;^ > 2 |l 6 l, it follows that ak,t = 2 {n-kL)\kl)'^ monotonically decreases with increasing 
L < Hence from ® it follows that lim — decreases with increasing L. ■ 

Lemma [T] shows the interesting result that the SNR gap between the residual CFO scenario 
and the zero CFO scenario decreases with increasing frequency-selectivity (L) of the channel. 
For the residual CFO scenario, with 7 ^ = we note that SINRfc[f] depends only on c and L, 
for sufficiently large M and H in ([7]) vanish with M —)■ 00 ). From O it is clear that 

for fixed c, the MSB decreases with increasing L and hence SINR^ [f] would increase. Similarly 
from d?]) we also note that for a fixed L and decreasing c, the MSB increases and therefore 
SINRfc[f] decreases. Hence, for a fixed desired information rate, i.e. fixed SINRfc[f], we must 
decrease c with increasing L, i.e., 7 ^ = must be decreased. Therefore the SNR gap with 
the zero CFO scenario decreases with increasing L (since with = 0, it is clear from (|7]) that 
SINRfc[f] is independent of L). This conclusion is also supported in Fig. |2l 


IV. Numerical Results and Discussions 


In this section, through Monte-Carlo simulations, we study the variation in the minimum 
required received SNR 7 ^ with increasing per-user information rate for different values of L (fixed 
N, K, M and iV„). We assume the following: carrier frequency fc = GHz, a maximum CFO 
of nfc (k = 0.1 PPM) and communication bandwidth = 1 MHz. Thus \ujk\ ^ 271 = 2 ^- 

At the start of every CFO estimation phase, the CFOs {k = 1,2,..., iC) assume new values 
(independent of the previous values) uniformly distributed in [— 2 ^) 2 Hm]‘ duration of uplink 
is Nu = 2000 channel uses and pilot length for CFO estimation N = N^. The PDP is the same 
for all UTs and is given by ah^i = 1/L, I = 0,1,..., L — 1; k = 1, 2,..., Ff. The information 
rate is computed using d?]), with = E[(cDfe — tokY] replaced by its expression in dU) with 
Gfc = 1 (see the line before Remark [T]). 

In Fig. [2] we depict the variation in the SNR gap between the residual CFO scenario and the 
ideal/zero CFO scenario, with increasing desired information rate Ik for the first UT {k = 1) for 
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Fig. 2 Plot of the SNR gap between the residual CFO and ideal/zero CFO scenarios with increasing information 
rate, Ik for the first user {k = 1), for different L = 1, 3, 5,10 and 20. Fixed parameters: M = 160, if' = 10 and 
N = Nu= 2000. 

M = 160 BS antennas and L = 1, 3, 5,10 and 20. Note that with L = 1 for = 3 bpeu, the SNR 

gap is ~ 4.22 dB, whieh quiekly deereases to a small value of ~ 0.07 dB when L = 20. This 

supports our eonelusion in Lemma [T] that the performanee of the TR-MRC reeeiver in massive 

MIMO systems with CFO estimation/eompensation proposed in [[6l, approaehes the zero CFO 

seenario performanee limit with inereasing frequeney-seleetivity of the ehannel. 
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